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Abstract 

We study the effect of quark and nuclear matter mixed phase on 
the radial oscillation modes of neutron stars. For this study we have 
considered recent models from two classes of equations of state namely 
the relativistic mean field theoretical models and models based on real- 
istic nucleon-nucleon interactions incorporating relativistic corrections 
and three body nuclear interactions. 
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1 Introduction 



Recently there have been important developments in the determination of 
neutron star masses in binary systems. Where as the best determination 
of neutron star masses, all in the range 1.35 ± 0.04 Mq (Thorsett 1999), is 
found for binary radio pulsars, limits on the masses of neutron stars from 
measurements of kilohertz quasiperiodic oscillation frequencies in low mass 
X-ray binaries lie in the region of 2Mq (Miller 1998a; 1998b; 1998c; see also 
van-der Klis 1998). In addition several X-ray binary masses have also been 
found to be high viz. Vela X-1 and Cygnus X-2 with masses (1.8 ± 0.2)^^© 
and (1.8 ± 0.4)Mq respectively (van Kerkwijk 1995; Orsoz 1999). These 
large masses put severe constraints on the equation of state(EOS) of dense 
matter and only very stiff EOS's are capable of sustaining such high masses. 
There are esentially two classes of EOS's for dense matter, one class based 
on relativistic mean field type of models with coupling constants treated as 
parameters to be fitted to observable quantities and the other based on re- 
alistic interactions between constituents obtained by fitting scattering data 
and using techniques of many body theory to evaluate energy and pressure. 
Both classes of models however, suffer from lack of sufficient experimental 
data. This situation is likely to be ameleorated by our understanding of the 
physics of heavy-ion collisions in laboratory experiments in near future. In 
the core of neutron stars where densities rise few times the normal nuclear 
densities, the state of matter is essentially unknown and the possibility of 
a phase transition to constituent quark matter or to hadronic matter con- 
taining hyperons and/or meson condensates has been extensively studied in 
the litrature (see for example Heiselberg 2001 and references therein). All 
these phase transitions whether first order or second result in softening of the 
EOS and thus aggravating the maximum mass limit. In a phase transition to 
quark matter, the neutron star may have a core of pure quark matter with 
a mantle of nuclear matter surrounding it and the two phases co-existing in 
a first order phase transition. Alternatively one may have a so called hy- 
brid star discussed by Glendenning (Glendenning 1992; ibid 1997; see also 
Heiselberg 2001), wherein the quark and nuclear matter coexist in mixed 
phase with continuous pressure and density variation - a situation obtained 
by applying Gibb's criteria to a two-component system. 

One important aspect of various compact stars is the study of their radial 
modes of oscillation. Radial oscillations give information about the stability 
of the stellar structure under consideration, and arc thus quite important 
in distinguishing between various models of the stellar structure. Also as 
is well known, radial oscillations do not couple to gravitational waves; as a 
result the equations governing radial oscillations are quite simple and it is 
relatively easy to solve the eigenvalue problem that leads to the descrete set 
of eigenfrequencies for the system. The eigenfrequencies form a complete set 
and hence it is possible to describe any arbitrary periodic radial disturbance 
as a superposition of its various eigenmodes. 
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For neutron stars, beginning with Chandrasekhar (1964), radial modes 
have been investigated in the hterature for more than thirty years and for 
various nuclear matter equations of state, (see for example Chanmugam 1977; 
Glass and Lindblom 1983; Vath and Chanmugam 1992; Kokkotas 2001) . In 
view of the fact that only a few EOS viz. those of Wiringa, Fiks and Fab- 
rocini (Wiringa 1988) using A14+UV11 and U14+UV11 interactions and 
the recent EOS of Akaml, Pandharipande and Ravenhall (Akaml 1998) in- 
corporating relativistic corrections and three- nucleon interactions in a new 
nucleon-nucleon interaction model (the Argonne A18 potential) arc prac- 
tically the only viable neutron star models capable of giving neutron star 
masses greater than 1.9 Mq, we have picked them for the present study. In 
addition we also consider the EOS given by Glendenning from the class of rel- 
ativistic mean field theoretical models (Glendenning 1997). These equations 
invariably allow a phase transition in the core. In this paper we study the 
effect of a mixed quark-nuclear matter core on radial oscillations of hybrid 
stars. 



2 Radial Oscillations of a non-rotating Star 

The equations governing the radial oscillations of a non-rotating star, using 
static, spherically symmetric metric 



-e'^df + e^^dr' + r\d0' + sin'Odct)') 



(1) 



were given by Chandrasekhar (Chandrasekhar 1964). The structure of the 
star in hydrostatic equilibrium is determined by the Tolman-Openheimer- 
Volkoff equations 

dp —G{p + p){'m + inr^p) 
dr r^fl — 



2GM 



47rr^p 



dm 

dr 

du 2GM(1 + ^) 
dr 



(2) 
(3) 



r 1 - 



2GM^ 



(4) 



where we have put c = 1. Assuming a radial displacement 5r with harmonic 
time dependence 5r^e''^^ and defining variables ^ = — and ^ = r'^e~^^ , the 
equation governing radial adiabatic oscillations is given by 



dr^ dr 



where 



P + p 



P + P 

dr r 
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(6) 
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A is related to the metric function through 



e- = (l-H^) (9) 
and 7 is the adiabatic index, related to the speed of sound through 

7 = (10) 
P dp 

Equation(5) is solved under the boundary conditions 

((r = 0) = 5p{r = i?) = (11) 

where 6p{r) is given by 



Equation (5) with the boundary condition (11) represent a Sturm-Liouville 
eigenvalue problem for u'^ with the well known result that the frequency 
spectrum is descrete. For u"^ > 0, lo is real and the solution is purely oscil- 
latory whereas for u>'^ < 0, u is imaginery resulting in exponentially growing 
unstable radial oscillations. Another important consequence is that if the 
fundamental radial mode ujq is stable, so are the rest of the radial modes. 
For neutron stars ujq becomes imaginery at central densities Pc > p^"**'^"' for 
which the star attains it's maximum mass. For pc = p^rttica/^ ^^^e fundamental 
frequency coq vanishes and becomes unstable for higher densities and the star 
is no longer stable. There also exists another unstable point at the lower 
end of the central density, namely there exists a minimum mass for a stable 
neutron star and the frequency of the fundamental mode at the minimum 
mass again goes to zero. 
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3 Results 



The only information required to obtain the structure of the star and eigen- 
values of radial oscillation modes is the knowledge of the EOS. For a given 
EOS, equations (2)-(5) are solved numerically by standard techniques under 
the boundary conditions (11). While numerically integrating the equations 
for each EOS we make sure that the eigenfrequency of the fundamental mode 
goes to zero at the maximum mass of the star. We have also checked that 
the frequency also vanishes at the minimum stable mass too. For the EOS, 
as discussed in the introduction, we use two relativistic mean field theoretic 
models taken from Glendenning (Glendenning 1997) and two potential mod- 
els incorporating relativistic corrections and three body interactions given by 
Akaml, Pandharipande and Ravenhall (Akaml 1998). Both class of models 
admit of mixed quark-nucleon phase in the core and correspond to matter as 
below : 

RFTIH : pure confined hadronic phase for nB < 0.26/m~^, mixed con- 
fined phase at intermediate densities and pure quark- matter at nB > 1.17 fm~ 
Nuclear properties correspond to K — 240 MeV and ^ = 0.78. 

RFT2H : confined hadronic phase for nB < 0.26/m^'^, deconfined phase 
for nB > 1.0/m-^ mixed phase inbetween. K = 300 MeV and ^ = 0.78. 

The models RFTl and RFT2 correspond to pure hadronic phase occuring 
in the core of neutron stars without the existence of any quark phase. 

PMTIH : Incorporating Argonne A18 potential along with three body 
interaction (A18-I-UIX model of Akaml et al. 1998) corresponds to pure 
confined hadronic phase for nB < 0.1/m~^, mixed phase at intermediate 
densities and pure phase for nB > 0.49/m~^. 

PMT2H : Incorporating relativistic boost correction Sv along with three 
body interaction UIX* with Argonne A18 potential model corresponds to 
pure confined hadronic phase for nB < 0.1fm~^, and pure phase for nB > 
0.74/m-3. 

The models PMTl and PMT2 refer to the above two models with pure 
confined hadronic phase. 

The quark matter in the above EOS models is described by the Bag Model 
with niu = nid = 0, nig = 150 MeV, the Bag constant B^ = 180 MeV and 
as = 0. 

To illustrate the effect of mixed phase on neutron star parameters, namely 
mass-radius relationship and on the frequency of radial oscillations, we have 
taken one model each from the above two categories considered, since the 
mass fraction contained in the crust of the star is a small fraction of the 



total star mass (< 2%), we have used the earher results (Baym 1971; Lorenz 
1993; Pethick 1995) of the EOS for matter densities < 0.1/m-^ In fig 1 
we have plotted the M-R curves and find that the effect of the existence of 
mixed quark-nuclear matter phase in the core of neutron stars is to reduce the 
maximum mass. The effect is more pronounced for the relativistic mean field 
theory model than for the potential model. Whereas for the RFT models also 
radius corresponding to maximum mass decreases, for the potential models 
it increases instead. For smaller mass the two curves RFT2 and RFT2H 
(PMT2 and PMT2H) come very close to each other; this is to be expected 
since these points correspond to low central densities by which time the effect 
of deconfinement becomes negligible. In fig2 we have plotted the frequency 
of the fundamental mode and the next mode as a function of central density 
for the pure neutron and hybrid stars in the two categories and find that 
the frequency exhibits oscillatory behaviour in the case of a neutron star 
with a mixed quark-nuclear matter core. In tables (la)- (4b), we provide the 
numerical results for the star parameters and radial mode frequencies of the 
models discussed in the text. In each table we list the central density pc, 
radius R, mass M, red shift Z — {1 — 2^^)-2 — 1 of the stellar model and 
the frequencies of the first two radial models. The model above the stability 
limit is marked by a star. 



4 Acknowledgements 

We thank J. V. Narlikar for providing hospitality at Inter-University Centre 
for Astronomy and Astrophysics, Pune. One of us A. G. thanks University 
Grants Commission, India for partial financial support. 



TABLE la: RFTl 



n y 1 n^"^ 
yc ^ J-U 


R 

IV 


M 

IV J. 


7 
zj 




til 1 1 




(Km) 






(KU7) 






1 n fi*^ 

xu.uo 




u.oo 


QR 

U. t/U 


fi 99* 


25.35 


10.95 


1.550 


0.31 


0.45 


6.32 


15.22 


12.09 


1.496 


0.257 


1.52 


6.26 


10.65 


12.78 


1.394 


0.22 


1.89 


6.16 


8.79 


13.07 


1.314 


0.19 


2.09 


6.14 


7.00 


13.35 


1.187 


0.17 


2.38 


6.05 


5.83 
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3.61 


13.70 


0.608 


0.07 


2.89 
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3.07 
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0.471 


0.05 


3.04 


3.18 



TABLE lb: RFTIH 
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Z 
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31.68 


10.20 


1.45 


0.31 


0.89 
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10.42 


1.45 


0.30 


0.29 


6.12 


18.50 


11.18 


1.41 


0.26 


1.84 


6.31 


15.26 


11.56 


1.36 


0.24 


2.04 


12.60 


11.68 


12.18 


1.24 


0.20 


2.02 


6.65 


8.97 


12.85 


1.10 


0.16 


1.80 


6.26 


6.69 


13.38 


0.98 


0.13 


1.90 


5.64 


4.58 


13.58 


0.84 


0.11 


2.84 


5.49 


4.15 


13.62 


0.74 


0.09 


2.81 


4.86 


3.61 


13.69 


0.61 


0.07 


2.89 


4.01 
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TABLE 2a: RFT2 
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7.04 
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2.64 


5.066 


3.62 


13.97 
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2.60 


4.143 


3.07 


14.17 


0.495 


0.057 


2.51 


3.146 


2.55 
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0.344 


0.036 


2.05 


2.394 



TABLE 2b: RFT2H 



PcXlO^^ 


R 


M 


Z 




1^1 


(gm/cc) 


(Km) 


{Me) 




(KHz) 


(KHz) 


32.07 


10.31 


1.430 


0.303 


0.776 


6.27* 


27.40 


10.68 


1.434 


0.289 


0.213 


6.13 


20.88 


11.40 


1.425 


0.261 


0.851 


5.57 


17.39 


11.68 


1.409 


0.248 


1.426 


5.83 


13.27 


12.20 


1.346 


0.219 


1.723 


6.16 


10.92 


12.66 


1.268 


0.193 


1.721 


6.14 


8.60 


13.22 


1.160 


0.163 


1.641 


5.72 


6.48 


13.68 


1.041 


0.136 


1.795 


5.09 


5.26 


13.83 


0.958 


0.122 


2.262 


5.00 


4.10 


13.89 


0.807 


0.099 


2.643 


5.07 


3.62 


13.97 


0.651 


0.077 


2.603 


4.14 


3.07 


14.17 


0.495 


0.056 


2.515 


3.15 


2.55 


14.81 


0.344 


0.036 


2.054 


2.40 
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TABLE 3a: PMTl 
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6.52 
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0.65 
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0.45 


0.06 
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0.35 
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TABLE 3b: PMTIH 
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31.40 
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0.46 
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0.38 
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7.61 
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12.05 
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0.34 


3.05 


7.70 


9.07 


12.15 
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0.29 


3.54 
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1.56 


0.27 


3.71 


8.17 


7.41 
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1.33 


0.22 


3.73 


7.81 


7.30 


12.22 


1.18 


0.18 


3.49 


7.10 


6.33 


12.23 


0.92 


0.13 


3.54 


6.73 


5.95 


12.26 


0.82 


0.12 


3.46 


6.31 


5.39 


12.32 


0.70 


0.10 


3.36 


5.61 



TABLE 4a: PMT2 
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TABLE 4b: PMT2H 
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Figure 1: Plot of mass in solar mass unit Vs radius in Km. 




Figure 2: Plot of frequency in KHz Vs central density in gm/cc 
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